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The rounding of the charge density wave depinning transition by thermal noise is 
examined. Hops by localized modes over small barriers trigger "avalanches" , resulting 
in a creep velocity much larger than that expected from comparing thermal energies 
with typical barriers. For a field equal to the T = depinning field, the creep 
velocity is predicted to have a power-law dependence on the temperature T; numerical 
computations confirm this result. The predicted order of magnitude of the thermal 
rounding of the depinning transition is consistent with rounding seen in experiment. 
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The model |T[ of an incommensurate charge density wave (CDW) as a deformable medium 
has been successful in explaining many experimental results both qualitatively and quanti- 
tatively [0]. In this model, the CDW is an elastic medium subject both to impurity pinning 
forces and an external drive force. At zero temperature, there is a sharp depinning transi- 
tion: when the applied electric field E exceeds a threshold value Ex, the CDW depins, and 
the sliding CDW carries an electric current Jcdw- Fisher P] has proposed that the behavior 
near the depinning transition can be understood as a type of dynamic critical phenomenon. 
In particular, the CDW current jcDW should behave as jcDW ~ (-E — Et)'^, with ( a critical 
exponent. Numerical work has confirmed aspects of this picture in finite-dimensional 
models, and values for the exponent ( have been determined 0,0. 

This same numerical work has shown that, in order to clearly see the dynamic critical 
behavior, it is necessary to examine fields within a few percent of the depinning field Et- 
It has proven difficult to obtain clear experimental results for applied fields this near to 
the depinning transition. A major cause of this difficulty is that the transition between 
the pinned and sliding state is not perfectly sharp ; there is an apparent rounding in 
the jcDw(-E') relationship. Besides affecting the measurement of currents for E near E^, 
the rounding makes it difficult to determine the threshold field E^ accurately. In order to 
make a careful comparison between experiment and theory, it is necessary to understand 
the rounding of the depinning transition. 

Various explanations for this rounding have been proposed. One possibility is that exper- 
imental samples are macroscopically inhomogeneous, with local regions depinning at distinct 
applied fields E. The bulk jcDw{E) relation would then be smeared. However, the sharp 
narrow band noise seen in many samples 0, indicates that the whole CDW is sliding at 
a uniform rate (though the narrow band noise is a finite size effect, its presence in finite 
experimental samples indicates uniform velocity). A somewhat related explanation is that 
provided by phase slip |p in macroscopically homogeneous samples. Coppersmith has shown 
that there exist exponentially rare regions with atypically weak pinning. Phase slip oc- 
curring on the border of such regions allows for an excess current below the bulk depinning 
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field Et- The main difficulty in attributing the rounding of the transition to this effect 



is the rarity of these atypical regions |T0|. In particular, the magnitude of this effect will 
decrease as the exponential of a negative power of the impurity concentration. It has also 
been proposed [0 that thermal noise may round the CDW depinning transition. 

In this paper, I examine the effects of finite temperature T on the CDW depinning tran- 
sition. Recent work |Tl| has examined thermal effects at small field, far from the transition. 



where a creep velocity oc exp[(£'T)~'^] is expected, for some exponent /i. I study here instead 
the thermal effects at fields E ^ Et, by examining the barriers that prevent the forward 
motion of Lee- Rice domains (regions of a size such that pinning and elastic forces are compa- 
rable and that act roughly as single degrees of freedom [Q). At first sight, thermal effects in 
CDW's might be expected to be very small, since the thermal energies are ~ 10^ - 10^ times 
smaller than the typical barrier energies estimated from the magnitude of the threshold field. 
This is consistent with the thermal creep for small fields, E <^ Et, being extremely small. 
However, at T = 0, for fields just below the depinning threshold, many of the barriers to 
forward motion are much smaller than the typical barriers [§,0- If the field is increased 
slightly, the smallest barriers vanish, and spatially localized instabilities are induced, leading 
to the "jumping" forward of individual Lee-Rice domains. Numerical calculations at T = 
for fields below threshold show that the jumping forward of a single domain may cause an 
"avalanche", in which a large region of the CDW slides forward, as shown in Refs. 
At fields just above threshold, there is no stable configuration, and these triggering jumps 
and resulting avalanches lead to a positive CDW velocity. The argument I develop here 
estimates the current induced by the inclusion of thermal noise, which causes "hops" with 
effects similar to those due to the "jumps" resulting from an increase in field. These "hops" 
result in avalanches, as shown in Fig. 1. For E = Et, this argument results in a prediction 
for the CDW velocity (directly proportional to the current jcDw) of the form 

v{T,E = Et) r^T'^'^ (1) 

where r is a non-universal exponent, with the values r = 3/2 for the usual Fukuyama-Lee- 



Rice model, and r = 2 for a "ratcheted-kick" model 0. The form of this scaling (with 
r = 3/2) is the same as that predicted by Fisher in mean field theory using a related, but 
distinct, argument 0. I find the power law behavior (|I]) to be consistent with numerical 
results in dimensions d = 2,3. For thermal energies on the order of 10"'* of the typical 
barrier height, these numerical calculations show a broadening of the transition on the order 
of 0.5%, in (i = 2, 3, in order of magnitude agreement with experimental observations. 

The Fukuyama-Lee-Rice model Hamiltonian on a cubic (or square) lattice in scaled 



units is given by ||4|, |14|J12 



(ij) i i 

where (pi is the CDW phase at spatial site i, h is the impurity pinning strength. Pi is a 
random pinning phase, uniformly chosen in the interval [0,27r), and V{(j) — {3) gives the 
shape of the pinning potential. The first term of Eq. (^ represents elastic interactions 
between nearest neighbor pairs of sites {ij), the second term models pinning by impurities, 
and the last term represents the CDW polarization energy. I consider here two potentials 
V: the usual, smooth potential V = cos{(f)i — (5i) and a sharp "ratcheted-kick" potential 
which is, effectively, a sawtooth function of cpi — f3i. Evolution of the 0j in time t is 
given by overdamped equations of motion dcpi/dt = —{dTi/dipi) + rii{t), with Langevin noise 
r]i{t), {rii{t)rii{t')) = ATE^Slt — t'), so that temperature T is measured in units of the barrier 
energy scale 2Et set by the threshold field. For T = 0, the velocity v = {d(j)i/dt)- is zero for 
\E\ < Et and is positive for E > Et, with f ~ /'» for small reduced field f = {E — Et)/Et- 
As there is no phase slip in this model, the temporally averaged velocity must be uniform 
throughout the system, and there is a single threshold field at T = 0. For \E\ < Et, 
the polarization is defined as P = {(pi). Numerical investigation of the sub-threshold state 
P JT^JT3[] shows that the linear polarizability of the CDW can be understood as the sum of 
contributions from spatially localized modes. As the applied field is increased, the modes 
with the smallest linear relaxation rates become unstable (via a saddle- node bifurcation), 
resulting in localized, irreversible jumps in the CDW phase. These jumps trigger avalanches 



over regions which become larger as threshold is approached; these jump and avalanche 
combinations are the dominant contribution to the nonlinear polarizability of the CDW . 

Consider the CDW configuration {(pi} at the threshold field, E = Et, and temperature 
T = 0. The localized modes of the CDW configuration can be considered as (linearly) 
independent degrees of freedom ipk {k being mode indices), each with size approximately that 
of a Lee- Rice domain. Each ipk sees an individual potential Wki'ipk)-, due to elastic, pinning, 
and drive forces, with Y,Wk{;ipk) = 'H. When the CDW is in a stationary configuration, 
each ipk is at a local minimum of Wk- Expanding Wk in powers of ipk-, with ■i/'fc = at the 
infiection point separating the local minimum from the barrier to forward motion, gives (see 
also Refs. ||,|): 

W{i^k) = Ia^jI + Bki^k - {Xk + 1){E- ET)i^k + 0{^l), (3) 

where A is approximately independent of k with a value determined by the strength and 
periodicity of the pinning potential and Bk and Xk are positive; for the softest modes k, 
Bk approaches zero, while Xk-, which incorporates the linear polarizability of neighboring 
domains, is of order one. The potential W{ipk) changes in response to an increase of the 
applied field. If E is increased by a small amount 5E, minima of the local potentials vanish 
where Bk is smaller than 6E/{xk + !)• As a result, some ipk jump forward and trigger an 
avalanche. By destabilizing a local mode, the increase of the field E leads to a current pulse. 
I now argue that, using Eq. (^, the effects of hops caused by a finite temperature can be 
directly related to the effects of the field induced destabilization of local modes. 

It is useful to note two numerical results from T = simulations of CDW's in finite 
dimensions: the apparent insensitivity of the exponent ( to the shape of the pinning potential 
and the relatively quick relaxation of the velocity in the moving state. In the rather different 
cases of the smooth pinning potential and the ratcheted-kick potential, ( is found to be 
universal to within numerical error, for d = 2, even though the mean-field exponents are 
distinct. The difference in the mean- field critical behavior of the two models results from 
the details of how a jump occurs. The apparent universality of ( therefore suggests that the 
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details of how a single domain jumps is unimportant in determining the critical behavior of 
the velocity in finite-dimensional models. This conclusion is supported by the observation 
that the velocity equilibrates well before the configuration reaches a periodic state. These 
numerical observations imply that the velocity is dominated by the avalanches which are 
triggered by the jumps, with the critical behavior of the velocity independent of jump 
dynamics and the shape of the sliding configuration. 

These results suggest that thermally-induced hopping over small barriers will have the 
same effect on the velocity as raising the applied field by an amount which causes the 
destabilization of the same soft modes, if the time scale for the hop is much smaller than that 
for the evolution of the resulting avalanche. For potentials of the form Eq. (0), the barrier 
height Afc behaves as ~ bI^"^. Since a finite temperature T will cause rapid tunneling 
over barriers of height ^ T, thermal noise triggers avalanches at a rate corresponding 
to an applied field E with {E - Et) ~ T^/^, assuming Xk of order one [|15|. It follows 



that, for small T at an applied field E = Et, the velocity v, due to avalanches triggered by 
thermal noise, scales as ~ T^/^, where the exponent r = 3/2. For the case where pinning 
potential V{(j) — f3) is given by the ratcheted- kick form, it can be similarly shown that 
r = 2. Thermal effects are therefore non-universal in the shape of the pinning potential (it 
is probable however that the thermal effects are universal if the noise is due to finite random 
kicks, rather than Langevin noise |^). The comparison of the effects of small reduced fields 
/ and temperature T also suggests the scahng form 

v{f,T) = T'^/-B{fT-'/-), (4) 

which is exactly the form proposed by Fisher for mean- field theory [Q], with B{u) behaving 
as u'^ for large positive u and decaying very rapidly as m ^ — oo. 

I have used numerical simulations for dimensions ci = 2, 3 to check the scaling behavior 
predicted by Eqs. The pinning strength h = {2.5)d is chosen so that the size of a 

Lee- Rice domain is approximately one lattice unit. The computations were performed on a 
Connection Machine CM-2, utilizing 16K processors for approximately 35h. 
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Fig. 2 shows velocity f as a function of dimensionless temperature T on a log-log scale, 
for fixed field E = Et and a smooth potential V. The straight lines shown are not fits 
to the data, but have slopes predicted by Eq. (|1]) with values for the dynamical exponent 
of C = 0.63, 0.85 in d = 2, 3, respectively (estimated error bars for ^ in d = 2, 3 are 
0.07; the value for in (i = 3 is from simulations on the ratcheted-kick model), and the 
predicted exponent r = 3/2. There is very good agreement between the data and Eq. (P, 
for velocities v ^ 0.1 (this velocity scale corresponds to the crossover between the high field 
and critical behavior in T = simulations f^.) If the slope is allowed to vary, a best fit 
yields C/t = 0.40 ± 0.04,0.60 ± 0.06 in d = 2,3, respectively. Similar calculations for the 
"ratcheted-kick" model in c? = 2, 3 also show power-law behavior consistent with ([^), for the 
above values of (, but with r = 2. It follows that the relationship between v and ksT/h is 
non- universal; this is not surprising given the derivation of Eq. ([^), which, when Langevin 
noise is used, depends on microscopic properties of the model. 

In order to examine the (non-universal) scaling function B of Eq. (^) for a particular 
case, velocities have also been calculated for fixed temperature and varying field. In Fig. 

(0) , scaled velocity t>T~^/^ is plotted as a function of scaled reduced field /T~^/^ for the 
smooth pinning potential in d = 2, for the above values of ( and r. The scaling form Eq. 

(1) describes the data well. If ( and r are allowed to vary so as to produce the best fit 
to a single curve, I find ( = 0.6 ±0.1 and r = 1.6 ± 0.2 (with subjective error bars). For 
dimensionless temperatures T of order 10~^, this plot indicates that the rounding of the 
transition occurs over a range of ~ 1% in reduced field, in d = 2. In general, the width in 
field E of the rounding of the v vs. E curve, will be some constant times EtT^^"^, independent 
of dimension. Using this criterion and the data of Fig. (0), the width of the rounding is 
found to be ~ 0.5% for T = 10"^ in = 3. 

The scaling predictions of Eqs. (|1|,§), which are supported by these numerical simulations, 
can be directly applied to experiment. The temperature scale set by the typical barrier is 
of the order Eb = 2pVflr-E'tAcdw, where Acdw is the wavelength of the CDW, which 
determines the periodicity of the pinning potential, p is the CDW charge density, and Vflr 
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is the Fukuyama- Lee- Rice domain volume. For typical values of the parameters in NbSea 
(AcDw = 4A, Et = 100 mV/cm, p = 2 x lO^C/m^, l/pLR = (O.l/im)^), at 125K, the 
dimensionless temperature is T = 3 x 10"^, on the order of the smallest temperatures 
examined in the simulations reported here, so that the rounding in this case should be on 
the order of 1%. Generally, the materials parameters, such as the elastic constant (which 
affects ^flr) and charge density, are strongly dependent on temperature, making it quite 
difficult to verify the scaling predictions by comparing with experiments where the physical 
temperature is varied. The most direct method for checking the behaviors in Eqs. (|l],^ 
may be to control the dimensionless temperature by varying the concentration of weak 
pinning impurities. For weak pinning, with impurity concentration n^, the dimensionless 
temperature T oc Ej^^ oc n^^^'^^^'''^^^^ 

I am pleased to acknowledge discussions with D. S. Fisher, P. B. Littlewood, J. Sethna, 
and C. Myers. The computational work was conducted using the resources of the Northeast 
Parallel Architectures Center (NPAC) at Syracuse University. 
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FIGURES 

FIG. 1. CDW velocity v vs. time t for a 128^ sample, with {E - Et)/Et = 0.98 and T = lO"^ 
in dimensionless units. Noise-triggered avalanches are visible above the local thermal fluctuations. 
The inset shows the avalanches which occur in the interval between the two arrows; dark regions 
show where the phase <pi advances by more than tt. In most of the volume of the avalanches, the 
phase advances by an amount close to 27r, while in the light regions, the phase changes by much 
less than 1. 

FIG. 2. Plot of CDW velocity u as a function of dimensionless temperature T at the threshold 
field Et, with smooth pinning potential. The key indicates the sample size and number of samples 
averaged over. The lines show slopes given hy (/t for exponent values given in the text for d = 2 
(dashed) and d = 3 (dotted) . 

FIG. 3. Thermal rounding of the depinning transition for d = 2 with smooth potential: 
symbols show scaled velocity vT~'^l'^ vs. scaled reduced field jT~^l'^ for two different temperatures, 
using r = 1.5 and C, = 0.63. The solid line shows power law behavior v ~ /^, which determines the 
scaling function B in the large fT~^/'^ limit. 



11 




Fig. 1 - A. A. Middleton, "Thermal Rounding ..." 
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